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Abstract 

Although the mirage mediation is one of the most plausible mediation mechanisms of supersym- 
metry breaking, it suffers from two crucial problems. One is the /i-/5/x-problem and the second 

£T) \ is the cosmological one. The former stems from the fact that the B parameter tends to be compa- 
rable with the gravitino mass, which is two order of magnitude larger than the other soft masses. 
The latter problem is caused by the decay of the modulus whose branching ratio into the gravitino 
^! pair is sizable. In this paper, we propose a model of mirage mediation, in which Peccei-Quinn 
Ph| symmetry is incorporated. In this axionic mirage mediation, it is shown that the Peccei-Quinn 
symmetry breaking scale is dynamically determined around 10 10 GeV to 10 12 GeV due to the 

l— ~~ '■ supersymmetry breaking effects, and the /z-problem can be solved naturally. Furthermore, in our 
model, the lightest supersymmetric particle (LSP) is the axino, that is the superpartner of the 

^ , axion. The overabundance of the LSPs due to decays of modulus/gravitino, which is the most 

CNJ ! serious cosmological difficulty in the mirage mediation, can be avoided if the axino is sufficiently 
light. The next-LSPs (NLSPs) produced by the gravitino decay eventually decay into the axino 
LSPs, yielding the dominant component of the axinos remaining today. It is shown that the axino 
with the mass of (9(100) MeV is naturally realized, which can constitute the dark matter of the 
Universe, with the free-streaming length of the order of 0.1 Mpc. The saxion, the real scalar 
component of the axion supermultiplet, can also be cosmologically harmless due to the dilution of 
the modulus decay. The lifetime of NLSP is relatively long, but much shorter than 1 sec, when 

^ 1 the big-bang nucleosynthesis commences. The decay of NLSP would provide intriguing collider 
c3 ■ signatures. 



1 Introduction 



Supersymmetry (SUSY) is one of the most promising candidates for new physics beyond the 
standard model (SM), which solves the naturalness problem associated with the hierarchy be- 
tween the electroweak scale and the Planck scale. However, since no super symmetric particles 
have been discovered yet, SUSY has to be broken. In order to prevent us from arising the new 
dangerous phenomenological effects, flavor-changing-neutral-current (FCNC) process and SUSY 
CP violation and so on, the SUSY-breaking mediation mechanism must be the one that does 
not give rise to these effects. In superstring theory, one of most plausible mediation mechanisms 
without these effects is the moduli mediation [I]. The contribution of the modulus, X, to the soft 
masses is proportional to its auxiliary component, Fx — m^, 2 /mx (in the Planck unit M-p\ = 1), 
where m 3 / 2 and nix are the gravitino mass and the modulus mass, respectively. On the other 
hand, in 4-dimensional supergravity (SUGRA), the mediation associated with the super- Weyl 
anomaly, the anomaly mediation [2] [I], also does exist, whose contribution to the soft masses is 
of the order of m 3 / 2 /87r 2 . Recently, KKLT [5] have proposed the interesting set-up to stabilize the 
modulus with the relatively heavy mass, nix — 47r 2 m 3 / 2 , by considering some non-perturbative 
effects. Therefore, in this and similar set-ups, the anomaly mediation contribution is comparable 
with the modulus mediation one. This type of mediation mechanism is often called the mirage 
mediation [6] [T5]. 

The mirage mediation is a natural mediation mechanism and quite interesting because it can 
solve the tachyonic slepton problem in the pure anomaly mediation, as well as has characteristic 
mass spectrums [2] [ID]. However, it has two crucial problems. One of them is the /i-problem. 
In the discussion of ref . [10] , although it is not difficult to obtain that \i is of the order of the soft 
masses, B becomes of the order of m 3 / 2 , in general, without fine-tuning of parameters. In the 
mirage mediation, since the gravitino mass should be O(10) TeV in order to obtain soft masses 
with the electroweak scale, this is problematic. 

The other problem results from cosmology. While in inflation the modulus field is shifted 
from the true minimum. When the Hubble parameter H decreased to the order of the modulus 
mass, mx, the modulus starts to oscillate around the true minimum with amplitude of order 
the Planck scale. The energy density of the coherent oscillation dominates that of the Universe 
before the modulus decay. If the modulus mass is around the electroweak scale, its decay produces 
immense amount of entropy after primordial nucleosynthesis and upsets the success of the big- 
bang nucleosynthesis (BBN). Such a cosmological catastrophe is known as the moduli problem 
[16]. One of the resolutions to the moduli problem is to invoke a relatively heavy modulus with 
mass of 10 5 GeV or larger. The heavy modulus then decays before the nucleosynthesis commences, 
which will not spoil it. Thus, KKLT set-up seems to be favored cosmologically. 

However, recently, it was found that non-thermal production of the gravitino from the modulus 
decay is sizable, which aggravates the moduli problem [T7] [33]. (See, for the case of the inflaton 
and the Polonyi field [IH] [22]). Such aggravation comes not only from abundance of the gravitino 
but also from the overclosure of the neutralinos if it is the lightest superparticle (LSP). The 
gravitinos with the electroweak scale mass decay into minimal supersymmetric standard model 
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(MSSM) particles after or during the primordial nucleosynthesis, and hence the decay products 
may spoil the success of BBN. This is called the gravitino problem. It is known that when the 
gravitino mass is as large as (9(10) TeV, there is no gravitino problem because of their sufficiently 
short lifetime. On the other hand, the neutralinos produced by the decay of the gravitinos are so 
abundant that the annihilation process is effective. However, its abundance after the annihilation 
exceeds the present dark matter (DM) abundance even if the neutralino LSP is the wino, whose 
annihilation cross section is most effective. Therefore, even if the modulus can decay before BBN, 
it causes the cosmological problem. This disaster may be remedied when LSP is lighter than the 
neutralino. 

We will solve these two problems simultaneously by the axionic extension of the mirage me- 
diation (axionic mirage mediation) . In SM, the instanton effect of the quantum chromodynamics 
(QCD) generate the following term in the Lagrangian 

£e = ^e e ^G^, (1) 

where g 3 is the strong coupling constant, G a the field strength of the gluon and 6 = 6 + 
arg(detA^) with the quark mass matrix M.. From the present experiment of the electric dipole 
moment of the neutron, 6 has to be smaller than 10~ 9 . In SM, however, there are no reasons 
that 6 should be so small, because it just a parameter of the Lagrangian. This is well known 
as the strong CP problem. One attractive solution to this problem is to introduce an anomalous 
global U(l) symmetry, called Peccei-Quinn (PQ) symmetry, U(1) P q [23]. When the U(1)pq is 
spontaneously broken, a pseudo Nambu-Goldstone boson, called axion a, appears. Since the axion 
couples with the gluon as 

£ = A (e + 4-) ^G^G* (2) 



64vr 2 V fa 

where f a is the decay constant of the axion, the potential of the axion is minimized at 6 + (a) / f a = 
0. Thus, PQ symmetry can solve the strong CP problem [23] ES] • From the above argument, it 
seems reasonable to consider the supersymmetric axion model, which is studied by many authors 
[27] [33]. In some supersymmetric axion models, the axino, which is a superpartner of the axion, 
can be LSP and become the DM candidate [21] [37]. It is realized by introducing only a Yukawa 
coupling between an axion superfield and messengers in the superpotential [39]. In such a model, 
the axion superfield is a flat direction at the tree level, and hence the axino mass is induced by 
quantum corrections. Thus, the axino can become lighter than the neutralinos. 

In this paper, we show that the axionic mirage mediation can solve not only the fi-/Bfi- 
problem but the cosmological moduli problem simultaneously. First, we briefly review the mirage 
mediation in section [2j In section [31 we discuss about our model and show that the PQ scale can 
be obtained within, so-called, the axion window by stabilizing the axion superfield. Throughout 
this paper, we will focus on the case where the axion superfield is stabilized at 10 10 GeV. The /i- 
problem is discussed in section HJ In this section, we also find that there is no SUSY CP problem 
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in our model. The cosmological implication is devoted to section [5] and O Since the modulus 
decay releases a huge amount of entropy, unwanted particles existed before is diluted away, and 
thus we discuss only after the modulus decay Although the modulus produces a number of 
gravitinos, it will cause no cosmological problem because the next LSPs (NLSP) produced by the 
gravitinos decay into axino LSPs. We consider the cases where NLSP is the bino, the higgsino, 
the stau, the stop or the wino, respectively. By discussing several processes to produce the axino 
by non-thermally, we find that the axino with mass (9(100) MeV produced such a way can explain 
the present DM abundance in any NLSP cases, if the annihilation process of NLSPs produced by 
the gravitino decay is not effective. The saxion, the real scalar component of the axion superfield, 
would be also harmless because the modulus decay dilutes its decay products since the lifetime 
of the saxion is shorter than that of the modulus due to the relatively small decay constant 
/pq ~ 10 10 GeV. We briefly comment on whether the axion can constitute DM when the axion 
superfield is stabilized at 10 12 GeV. Section [7] is devoted to summary. 1 



2 Mirage mediation of supersymmetry breaking 

First we briefly review the mirage mediation of the SUSY breaking [6] [H]. Let us consider the 
following supergravity (SUGRA) Q function and superpotential 2 : 

n = -3(x + xt) + (x + xt)^|g / | 2 , (3) 

W = W (X) + 1\ukQiQjQk, (4) 
o 

where X and Qj represent the modulus field and the chiral matter fields, respectively. The 
modulus superpotential, Wo(X), which is responsible for stabilization of moduli, arises from 
non-perturbative effects like gaugino condensation or stringy instanton as well as with a mod- 
ulus independent term from various origin. In the KKLT construction [5], it is given by Wo = 
Aexp(-bX) + C 3 . The function Q has the relation to the Kahler potential, K, via 

n = -Ze~ K '\ (5) 

From the SUGRA fl function and the superpotential, the chiral/gauge fields sector of the SUGRA 
Lagrangian can be written by using the compensator, $, 



C = I d 4 9 &$tt+ I d9 



+ ^f a (X)W aa W a a 



+ h.c, (6) 



where f a (X) is the gauge kinetic function. Assuming the shift symmetry, X — > X + iu, in absence 
of Wq(X), its form is restricted to f a (X) = p a X + q a . After integrating out F$, the scalar part 



introduction of a singlet field to solve the /j,-/£?/i-problem and the cosmological problem was also considered 
in [38]. The latter problem was solved by thermal inflation. 

2 Throughout this paper, we apply the Planck unit Mp\ ~ 2.4 x 10 18 GeV = 1 unless we explicitly mention. 
3 Here we consider the single modulus case only. Extention to the multi-modulus case is straightforward. 
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recovers the conventional SUGRA potential in the Einstein frame by choosing ($) = (e K ^), 



VsUGRA 



e K (K IJ D I WD J W - 3\W\ 2 ), 



(7) 



with iV($) = K^ 1 /3 + m 3/2 . Here, a SUSY configuration, F 1 = -e K / 2 K IJ DjW = 0, gives 
the extremum of Vsugra if it exists. It is often the case that this corresponds to the minimum of 
the potential and ends up with the AdS vacuum, Vsugra = — 3e x |IV"| 2 = — 3|m 3 / 2 | 2 . To achieve 
observed approximately Minkowski vacuum with positive cosmo logical constant, we introduce 
the uplifting potential, V = Vsugra + Vjif t ~ + . This additional SUSY breaking potential could 
be remnant of stringy effects like the anti-D3 brane in the original KKLT construction [S] or 
field theoretic "hidden sector" effects separated from the first term of Vsugra (F-term uplifting) 
[lQ]-[43j. Including such a modification, the initial SUSY configuration of the modulus field, 
DxWo(X) = 0, shifts by \5X\ ~ M Pi (m^ 2 /m x ) and this induces modulus SUSY breaking, 
\Fx/2Xr\ ~ m 2 , 2 /mx where Xr = Re(X). It is known that the non-perturbative stabilization 
predicts significant enhancement of mx/ m 3/2 due to derivative of the exponential factors in W ~ 
777.3/2- Thus Fx/2Xr is considerably suppressed relative to 777.3/2. This means that, in the visible 
sector, the anomaly mediated SUSY breaking of order (F$/(<3>))/87r 2 ps m 3 / 2 /87r 2 is equally 
important as the modulus mediated contribution. The relative phase of two sources of the SUSY 
breaking potentially causes the sever SUSY CP problem [H]. However, if the complex phases in 
Wq can be rotated away by symmetries which are explicitly broken by Wo {e.g. the R-symmetry 
and the shift symmetry), the relative phase in Fx / '(F$ / ($)) vanishes, which is the case in the 
KKLT set-up jl5] [9] [TO]. This new class of the SUSY breaking is dubbed "mirage mediation" 
after a peculiar renormalization group (RG) behavior of its mass spectrum which we will discuss 
later. 

Any direct coupling with the source of the SUSY breaking in the Kahler metric undermines 
such a loop suppressed effect in the scalar sector. Thus the sequestering of the SM fields and the 
hidden sector is essential for the mirage mediation. Actually, it is known that naive geometrical 
separation only does not guarantee the sequestering in string theory [17] [48] . However, strongly 
warped geometry, which is the heart of the original KKLT proposal, fulfills this requirement [5| 



In literature the following constant is defined to parametrize the relative strength of the two 
SUSY breaking contributions, 



The KKLT set-up with the uplifting potential generated by the anti-D3 brane on the tip of the 
warped throat predicts a = 1 + (9(l/47r 2 ), while there are variety of proposals which lead to 
different values of a = 0(1) ©-[ID] ED] S3- 

The soft SUSY breaking terms of the visible fields are parametrized as, 



[8] Hi] [49] [50]. 



a = 



7773/2 777 3 / 2 




(F X /2X R ) ln(M Pi /m 3 / 2 ) Aii\F x /2X R ) ' 




777/10/1 + A IJK Yij K ^j^j^ K + h.C, 



(9) 
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where A a denotes the gaugino and (pi is the scalar component of Qi in canonical normalization, 
while Yuk = ^ijk/ \f{X + X^) qi+qj+qK represents the canonical Yukawa coupling constant. We 
perform the redefinition of the matter fields, Qi = $Q/ to identify the effect of the anomaly 
mediation, 



d A 6 



d 2 6 



-3(x + xt)|$| 2 + (x + xt)^|g / | 

$ 3 W(X,Qj/$) + l -f a {X)W aa W a a 



+ h.c. 



(10) 



It is straightforward to derive the canonical gaugino mass from the above gauge kinetic function, 

Fx 



m Xa =F x d x \nRe(f a (X))=c Xa 



2X 
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while the other soft SUSY breaking terms are read off by substituting the equation of motion for 
the F-component of the light fields, Qj, into the above Lagrangian, 

2 „ 



m, 



Qi 



x 



2X 



R 



Mjk = (qi + qj + qx] 



'X 



2X 



(12) 



Here we used the fact that Xijk is modulus independent due to the shift symmetry. It is noted 
that this tree-level expression derived at the cut-off scale A of the SUGRA description does not 
have contribution from F<$> because of the classical invariance of the Lagrangian under the super- 
Weyl transformation. However, since the super- Weyl transformation is anomalous [2]-[3], we have 
F$ dependent terms at one-loop level via, 



fa(X) - f a (g, X) , (X + X^ - {X + X^Zj (JL x + x^j 



(13) 



where \x represents the renormalization scale. This is the anomaly mediation [2] [3] which is loop 
suppressed and works whenever F$ ^ even without the modulus contribution. Including this 
effect, it is noticed that the gaugino mass reproduces the original tree- level value at, so-called, 
mirage messenger scale, M mms = A(m3/2/Mpi) a / 2cA », which is independent of the choice of gauge 
group if c\ a is universal pU] . This is a consequence of cancellation between the quantum correction 
to the modulus mediation and the anomaly mediation. The scalar mass parameters also show 
similar behavior at one-loop level if the relevant Yukawa couplings are negligible or c\ a = 1, 
Qi + Qj + Qk = 1 are satisfied. The name of the mirage mediation is taken from the fact that the 
scale, M mms , does not correspond to any physical threshold, although the soft parameters show 
the systematic pattern (eq. fllll) . (fT2"|) ) at this scale, suggesting enhanced symmetry. 



3 The model 



Let us consider the hadronic axion model [25] in the above set-up by introducing the axion 
superfield S and N pairs of messenger fields \& and 1 l r , Qj = (S,^/,^/,Qi), where we use a small 
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subscript, i, to denote matter fields in MSSM. Here, S is completely singlet under any unbroken 
gauge symmetry and \l/ and ^ are vectorlike representation of the SU(5) gauge group. We assigned 
the PQ charge as Qpq(S) = —2, Qpq(^) = Qpq(^) = 1 and Q P q(X) = 0. The superpotential 
has Yukawa coupling, 

w = \sm + ---, (14) 

which is allowed by the PQ symmetry. S is a flat direction in SUSY limit, however, it is lifted by 
the SUSY breaking, Fx,$, caused by the mirage mediation. We assume that S is stabilized far 
away from the origin, breaking the PQ symmetry. This gives a huge mass for the messengers, so 
that they are integrated out at \i = \(S)\ ~ /pq, leaving S dependence in / and Zj at low energy 




(15) 



This introduces interactions among saxion/axino and the SM fields suppressed by /pq, which is 
much stronger than those by Mp\. Later we will see these interactions play important roles in 
cosmology. 

The scalar potential of S is derived from the Kahler potential, 



C = J d*6(X + X*) k Z s U^X + xA\S\ 



(16) 



where Z$ is the wave function renormalization of S at fi = \S\ and k = q$. We can obtain the 
relation among F-terms from the equation of motion of F$ as 



Ei ~ k Fx 

S 2X R 



2^1n|S'| 



(17) 



Here and in what follows, we use an abbrebiation, / ($) — ► F$ for the sake of brevity. It is 
noted that Fx/2Xr is smaller than F$ by one-loop factor in the ballpark, while din Zs/ din \ S\ = 
7s/(87r 2 ) corresponds to the anomalous dimension of 5* at [A — \S\ suppressed by the same factor. 
Integrating out the auxiliary component of S, the scalar potential is obtained as 



v(\s\) 





F x 




2X R 



1 d 2 lnZ s 
4d(ln\S\)< 



\F*\ + 



d 2 InZ, 



2 \dXdln\S\ 



F*tF x + h.c. ^|^| 2 = m||^| 2 . (18) 



The extremum of this potential is determined by the condition, 



d\S\ 



d 

~dnV\ 



rrii 



\S\=0. 



(19) 
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If m 2 s (fi) crosses zero at some scale \x = //q, we expect this condition holds in the vicinity of 
\S\ = Hq since dm 2 s /d\n \S\ ~ m 2 s /8ii 2 <C m|. In the case of the deflected anomaly mediation [3S] 
which does not include the modulus contribution we have introduced, the S potential is given by 
the second term in the parenthesis of eq. f|T8|) [30J, 

" s|S |Ft|2 ** ~ ( is?) ' N [ms) " 5{5N + 2)A2(S)] m ]Ftf ' (20) 

where we have neglected the gauge coupling constants except the strong coupling, g 3 , and A is 
normalized canonically as X 2 (S) = X 2 /(X + X<) k+q9+q * ZsZyZ^, . It is known that this form does 
not change by the RG evolution 4 . At a large value of S, the A 4 term dominates and it leads 
to m| > 0, while at a small value, m 2 s seems to become negative because of the asymptotic 
freedom of #3. Therefore at first sight S is expected to be stabilized at some scale {S) between 
the cut-off and the electroweak scale. However, in the case of the deflected anomaly mediation, we 
needed more than 3 flavors of messengers in order to avoid the tachyonic slepton [39] . Since the 
introduction of many number of messengers spoils the property of the asymptotic freedom, such a 
stabilization mechanism will not work without introducing new gauge interactions other than the 
SM. On the other hand, in our case, it is enough to introduce only one or two messengers because 
the tachyonic slepton problem has already been solved by the modulus mediation contribution. 
Furthermore, the positive contributions to the messenger masses, m 2 ^ ^ and A-term, A s ^ from 
the modulus mediation (eq. (fT2l) ) drive m 2 s to negative more strongly than those of the deflected 
anomaly mediation as read from the RG equation, 




Thus, the potential eq. fllSp may have a minimum at a scale (S). Actually, if c\ a = 1, k = and 
+ Qqi = 1 ar e satisfied, m 2 s crosses zero at M mms irrelevant to the other parameters in the model 
as discussed in the previous section [10] . According to the RG equation this is at least a local 
minimum as far as m% + m% > 0, which is likely satisfied with the enhancement by the strong 
interaction. If we increase (or decrease) + q^, breaking the mirage condition, the minimum 
will sit above (or below) M mms . In fig {1] we show the S potential for the gauge kinetic function, 
/ oc X (c\ a = 1) with N — 1, k = 0, q% — = £ — 1/2 and a — 1. The minimum comes slightly 
below M mms due to the one loop correction in the derivative of the S potential. 

The PQ scale /pq ~ (S) is constrained by various astrophysical and cosmo logical considera- 
tions. The relic abundance of the axion is given as [52] 

a _f^ (lo^vf (!) 2 ' < 22 > 

4 Precise treatment requires a separation of the triplet and doublet components in ty/if? and A. We take this 
into account in the following numerical calculation. 
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where G is the misalignment angle and h the Hubble constant in unit of 100 Km/sec/Mpc. Apart 
from our ignorance of 6, to avoid the axion overclosure, /pq < 10 12-13 GeV. On the other hand, 
the cooling of the SN 1987A puts a lower bound, /pq > 10 9 GeV [53]. From these observations, 
/pq should be in the range of 10 9 GeV < /pq < 10 12-13 GeV. Let us examine whether this bound 
can be satisfied by our model or not. In figfS], we depict (S) as a function of a for c\ a — 1, N — 1, 
k = and three rational choices of I = 0, 1/2, 1. In the KKLT set-up with the uplifting by the 
anti-D3 brane (a ~ 1), the SN bound on /pq is marginally satisfied with I = 1/2 because m| 
crosses zero at M mms ~ 10 10 GeV for a ~ 1. While stepping up to £ = 1 the observed DM could 
be saturated by axion taking into account involved ambiguities. It is plausible that i sits in the 
range < t < 1 in the string effective theory because the real part of the modulus represents the 
volume of the extra-dimension where the gauge fields propagate and i is the scaling dimension 
against this volume for the matter kinetic term which couples to these gauge fields. Combining 
this with the axion window we obtain a constraint, 0.3 < a < 1.5 in this set-up. We numerically 
checked that these results are robust against changes in N and A because both m|(A) and its 
RG equation depend them mostly through the overall factor iVA 2 . Once we turn on k > 0, the 
modulus contribution to m|(A) becomes competitive and defers its crossing zero in RG evolution, 
leading to lower / P q. Thus k = is only allowed choice for a ~ 1 unless A is sufficently large (> 1) 
to dilute the modulus contribution. To push up /pq for k > or a > 1.5 (e.g. [U] [53] 
we may need to introduce new gauge interactions which couples to the messengers and enhance 
m %/q, significantly via RG evolution. 

Next we estimate the physical masses of the component fields in 5*. The scalar component is 
decomposed into the saxion a and the axion a, so that 5* = /pQexp[(cr + ia)/v^2/pQ]. The mass 
of the saxion a at the minimum of the potential is given by, 



, ld>V(\S\) 



2 d\S\< 



S={S) 

dm 2 s 5NX 2 




2 + mi 



dln\S\ 8tt 



2 



dln\S\J \ dln\S\ 
(ml + m\ + \A X \ 2 ) , (23) 



where we used eq. ( 1191) to derive the second line. Assuming that m|, and ml are in the same 
order of magnitude with other superparticles, m 2 is one-loop suppressed against the other soft 
breaking mass squared. From eq. f)16p . the mass of the axino a, which is the superpartner of the 
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axion, is calculated by 



C 



J d A 6{X + X r ) k Z s 



J d i 6{x+xmz s + l - 




,1 + IM |S| 



, &S 1 
In ^-r^ H 



d 2 Z, 



1 <9 2 lnZc 
2dln|S|3Xt 



5| 
5 



a In 1 51 
1 <9 2 lnZ 5 



<9(ln|S|) 2 
5t 



111 



sts 
$t$ 



9Zs 

dX 



F x 9 2 + h.c. 



a a, 



4«9(ln|S|) 2 * \ ,S' 

where in the last line we have normalized the fields canonically. Thus, the axino mass is 

1 



rria 



8tt 2 
5N 



2 9Xt 



4-^(|.s*i)/-: 



s 



5t 



(24) 



(25) 



where 75 = (djs/dlnfi) and Ag^f is the sum of the A terms of the mirage and the anomaly 
mediation. One can find that the axino mass arises at two-loop order in contrast with the masses 
of the other fields being at one-loop. Thus it will be LSP and candidate for DM. We discuss later 
how heavy is the axino which can explain the present DM abundance. 



4 The \i- 1 B /i-problem 

An important property of this model is that it provides a natural solution to the /i-/5/i-problem 
(Here, we used the same notation \i for both the higgsino mass and the renormalization scale). To 
obtain the phenomenologically viable model, the //-term and the S/x-term should be generated 
at the soft mass order. Let us consider the following superpotential and fl function (See, for 
instance, [30] [39] ) 

W = y 1 TH 1 H 2 + y 2 S 1 S 2 T, (26) 
n = |Si| 2 + \S 2 \ 2 + \T\ 2 + kS\S 2 + h.c, (27) 

where we introduced new singlets Si, S 2 and T, whose PQ charge is assigned as Qpq(Si) = 
Qpq(S 2 ) = —2 and Qpq(T) = +4, and y±, y 2 and k are constants 5 . For simplicity, the modular 
weight of singlets is set to be zero. The assignment of PQ charge to new singlets leads to 
Qpq(Hi) = Qpq{H 2 ) = —2. This implies that we cannot obtain the /i-term from the following 
Kahler potential: AK = £(X + X^)~ n HiH 2 with constant £ [57J. The last two terms in eq. (!27|) 

5 Here, one can interpret the new singlet 6*1 or S2 as the previous one S. 
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are, in general, cannot be forbidden by the PQ symmetry 6 . As was discussed in the previous 
section, the Yukawa coupling of the messengers with Si yields a non- vanishing VEV for Si. Then, 
S*2 and T become massive and we can integrate out T. Its equation of motion leads to 

S 2 ^ -yi^k. (28) 
Substituting eq. (j2S|) into eq. f[2"Tl) . we can obtain 

/ct 
d A 9 -±HiH 2 + h.c. (29) 
Si 

Thus, the operator eq. fl29|) generates the /i-/-B/i-term in canonical normaization: 



F ] 

B = « yiiA 



Si , („ , „ \ r * 



-^ + (qh 1 +te 2/ . )V 

. Si 



(31) 



where we neglected field dependence in k and sub-leading terms from Z Hl ^ H , 2 7 . Since they are 
the same order of soft masses, the ^-/i^/i-problem can be solved. It is noted that in our model 
there is also no SUSY CP problem. Since the phases of two SUSY-breaking F-terms, F$ and F x , 
can be aligned [8] [9] [33], we find that the phase of B, which is obtained as 

Fa Fx Fx Id In 

B = —£- (to + ^(k~ q Hl ~ q H2 )^- ~ ? — (32) 

Si **r ^am|S| 

can be rotated away simultaneously with that of the gaugino mass and A-term. Therefore, we 
constructed a model in which there are no /^-problem as well as SUSY CP problem, and it leads to 
the axino LSP which is lighter than the neutralino. The appearance of the lighter LSP than the 
neutralino is favorable for the cosmological point of view. As we mentioned earlier, the moduli 
problem may be solved if there is the LSP which is lighter than the neutralino. 

So far, we have assumed that the singlet is stabilized by messengers as explained in the previous 
section. However, now the extra singlets can play the role of messengers, and the light quarks 
have PQ charges [26]. Thus, with sufficiently large Yukawa couplings, we can also construct a 
model without messengers. 



6 When we write down in the superpotential all terms that are allowed by the PQ symmetry with minimal 
Kahler potential, we find that those terms can be translated into the form of eq. and eq. ([27)1 . We also find 
that the coefficient k does not necessarily become small when y\ and yi come to be small. 

7 This corresponds to the radiative correction to n and B/j,. Note that if yi t 2 is sufficiently small the above 
results correctly give the input of (virtual) RG running at the unification scale in MSSM, although these terms 
are actually generated at /x = (Si). 
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5 Cosmology 



In this section, we discuss the cosmological implications of our model, in particular we shall 
closely investigate how it solves the LSP over-abundance problem due to gravitino decays. Since 
contributions to SUSY breaking from S and X are suppressed, we need to introduce an additional 
source of SUSY breaking, in order to cancel the cosmological constant. In KKLT, SUSY is broken 
by an anti-D3 brane. Otherwise an additional field Z is introduced. If supersymmetry breaking 
is of dynamical origin, its mass will be much higher than the electroweak scale. In what follows, 
we assume this is the case. During the inflationary epoch, the scalar fields such as Z, X and 
a are deviated from their true minimum with some magnitudes. After inflation, each of fields 
will start a coherent oscillation when the Hubble parameter, H, becomes comparable to its mass 
except for the saxion. The behavior of the saxion is exceptional, which we will discuss shortly. The 
oscillation energy behaves like a matter, whose energy density is red shifted more slowly than that 
of radiation, so that it may dominate the energy of the Universe if the initial amplitude is not so 
small compared with the Planck scale. It is plausible that the modulus has the initial amplitude 
of the order of the Planck scale, while that of Z is not the order of the Planck scale but the 
intermediate scale much smaller. Through out this paper, we assume that Z decays sufficiently 
fast not to dominate the Universe. The coherent oscillation of the modulus would commence 
before the reheating due to the inflaton, <fri, decay is completed, if the reheating temperature is 
not so high. On the other hand, the time when the saxion starts to oscillate depends on whether 
the saxion is trapped at the origin or deviates from the origin during the inflation. 

First, we consider the case where the saxion field is trapped at the origin during the infla- 
tion. This happens if the Hubble induced mass squared for the saxion at the inflationary era is 
positive. After the inflation, the Universe is reheated. The thermal effect through the messenger 
fields generates the temperature dependent effective potential for the saxion field. The resulting 
positive mass squared proportional to the temperature squared confines it at the origin until the 
temperature T goes down to T = Tc ~ 1 TeV where the thermal effect ceases to dominate over 
the zero-temperature contribution eq. fT23|) . The ratio of the energy density of the saxion to the 
entropy density at this time is 

p_a_ ^ 45 ml(S) 2 

s ~ 4tt 2 g^Tl 

~ 1 1 x 10 13 GeV l^EAY (J<L-Y ( m ° V ( ffi V (33) 

_ 1.1 x 10 OeV^ ioQ j [ 1Tey ) UoGeW \lG* GeV J ' { } 

On the other hand, since the modulus commences the coherent oscillation at the time when 
H ~ rrix, 

Px 45 m\Ml x 



s 4vr 2 ^(m x M P1 ) 3 / 2 

« -^GeV^H^f. ,34) 
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The ratios remain constant until these fields decay unless other entropy productions occur. Com- 
paring eq.(J2SD with eq. (1MI) . we find the energy density of the saxion is much larger than that of 
the modulus, and hence the saxion dominates the Universe. The saxion mainly decays into a pair 
of axions with the decay width 

1 771^ 

T CT ~ I> -> aa) = A 7 . (35) 

Since the axions produced by the saxion behave like neutrinos, they increase the Hubble expansion 
rate, and hence, the abundance of 4 He. Thus, once the saxion dominates the Universe it would 
upset BBN unless there is an extra entropy production {e.g. thermal inflation [58]) after its decay. 

Next, we consider the case where the saxion is displaced from the origin during the inflation. 
In this case, the saxion has at most the initial amplitude of the order of Mpi, S in < Mpi, and 
will start to oscillate when H ~ m a . We should note that the thermal effect generates the 
effective potential with positive curvature around the origin. The effect becomes irrelevant for 
the field value of saxion larger than the temperature because the messenger fields acquires masses 
larger than the temperature and the contribution to the effective potential is suppressed. How 
the saxion zero-mode moves under this thermal potential would be complicated, investigation of 
which is beyond the scope of this paper. Being aware of the possibility that the saxion field may 
be trapped at the origin, we assume that this is not the case and the saxion field follows simple 
damped-coherent oscillation with the initial amplitude of the order of the Planck scale. After the 
reheating of (pi, since the oscillation energies of X and a soon dominate the energy density of the 
Universe because of its relatively long lifetime. Their energy densities after the reheating are 

P^)^I T / X 2 P^l^^Sl (36) 
s 8 R m s 8 R in v ; 

where t\ is the time when the inflaton decays and the reheating temperature of the inflation. 
Here we consider the case both X and a start coherent oscillation during the era of inflation 
coherent oscillation (before the reheating of the 0/). However, since the decay constant of the 
saxion is much smaller than the Planck scale, the saxion decays much faster than the modulus. 
In fact, the lifetime of the saxion and the modulus is (See, eq.( |35|) and eq.(|52l)) 

"V y 3 ( 0) 



^ 13x10 sec WcSvJ (iraj< (37 > 

^OxlO-secy^)" 3 , (38) 

where we have set d g = 1. Therefore, even if the saxion has initial amplitude of the order of 
the Planck scale, the saxion decays faster than the modulus, and hence the modulus dominates 
the Universe as long as the confinement of the saxion in the thermal potential does not occur. 
When the modulus dominates the Universe, its decay produces a large amount of entropy. Let 
us estimate how large entropy is released by the modulus decay. When we denote that the 
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time at which the coherent oscillation of the modulus commences, the ratio of the energy density 
of X and that of the radiation at the modulus decay time, t x , is 



PR(fi) Pnwrn) p*(£c)\t?j -6 Xta \T?) ' (39) 

where T x is the decay temperature of X and R the scale factor. Therefore, unless the initial 
amplitude of X is not much smaller than the Planck scale the entropy increase factor A is 

when, for instance, ~ 10 6 GeV and Tf ~ 0.1 GeV 8 . This relatively large entropy increase can 
sufficiently erase unwanted particles produced before the modulus decay. The axions produced 
by the saxion decay are also diluted. The ratio of the energy density of the axion, p a , to that of 
the radiation at the modulus decay can be estimated as 

eMI _ Paw m) „ (iA 2/3 „ 3 x 10 - 3 (41) 

where t a d is the decay time of the saxion and we have used eq. (!37I) and eq. (l38l) in the last step. 
In order that the produced axions do not change the Hubble expansion rate, its energy density 
should be less than that of the one neutrino species: 



Po_ 

PR 



< —. (42) 
43 

lMeV ^° 



Comparing eq. (f4T!) and eq. (T42l . the axions produced by the saxion would not spoil BBN. There- 
fore, we find that it is sufficient to discuss the cosmic evolution only after the modulus decay in 
the case where both the modulus and the saxion have the initial amplitude of the order of the 
Planck scale. In what follows, we address the implications of the modulus decay to cosmology. 

The modulus mainly decays into gauge multiplets through the gauge kinetic function. (See, 
the second term in the last line in eq.flHJ)-) The decay width into gauge boson pairs and into 
gaugino pairs is |T7] [18] 

r( X^ 9Hr( ^ AAH J_(^) 4 2|, (43 ) 

where Nq is the number of the gauge bosons, which is 12 for MSSM. The numerical factor of 
order unity, d g , is defined by 



-i 



dX / 



(44) 



3 The modulus decay temperature with mx — 10 6 GeV is of the order of 0.1 GeV. See below. 
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where the subscript X represents derivative with respect to X and the real part = Re/, 
respectively. The function G is the total Kahler potential defined by G = if+ln IW^ 2 . Remarkably 
it was recently recognized that the decay of X into the gravitino pair is not suppressed. The decay 
width is given by [17] [18] 9 

r ( x^3 /2 <M = j£g, (45) 

where d 3 / 2 is also order one coefficient defined by 

2 

777 

(Gxxt)- 1/2 (e G / 2 G x )= C / 3 / 2 ^. (46) 

The modulus also decays into the saxion pair and the axino pair. The relevant interaction for the 
decay into the saxion pair is 



C-X<T(T — — - 



k dFx^ 1 d 2 \nZ s dF% 
(2X R f X dX 2dX^d\n\S\ * dX 



SXaa, (47) 



where 5X = X — (X). In eq. (1471) . the first term gives the dominant contribution because 

J 3/2/ 



{dFx I dX) ~ mx and (F x ) — mi /2 /mx- The decay width can be estimated as 



v ; 128tt \m x J M|, V 



The interaction of X with axino comes from eq. (j2 



1 d 2 \nZ s ^F x t 

£xhh ~ ——bXaa. 49) 

2<91n|S|«9Xt dX 1 ; 



Since the derivative of the anomalous dimension with respect to X is 

1 d 2 lnZ« 5NX 2 



2 01n|S|0Xt 16tt 2 ' 
the decay width into the axino pair is 

5NX 2 \ 2 m 



(50) 



IYX a a) ~ — ^ ( ) ^4. (51) 

; 128vt V 16tt 2 J M 2 X y 1 



9 In some cases where Z is light and has a minimal coupling, the mixing between Z and X might play an 
important role [59] [60], suppressing the decay of X into the gravitinos. However since mz 3> mx, this does not 
work in our case. 
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Comparing these partial decay widths, we find that the modulus dominantly decays into gauge 
boson/gaugino pairs with the total decay width 



Ng \ ,2 m x 

X ~ 16tt I 12 ) 9 M? ' 



S ) ■ (52) 

pi 



Thus, the reheating temperature of X is 

90 ^ 1/4 



where g*{T^) is the effective degrees of freedom of the radiation at the reheating. Before discussing 
how the present DM abundance is composed, we should comment on saxions produced by the 
modulus decay. The branching ratio of X into the saxion pair is obtained from eq. (j48l) and eq. (!52|) 



The yield of the saxion produced by the modulus decay is 

yX £ ±_R_ nX 

° ' 2 m x a 

- ^ U d g \ 10 J V 12 J Vl06GeVy llO^GeVy- 

However, from eq.( !37l) . we find that the saxions decay immediately when they are produced by 
the modulus, and hence their decay does not spoil the success of BBN. 

In our model, as we have already mentioned, the axino is the LSP and the DM candidate. 
Thus, in the following discussions, let us estimate the axino abundance and discuss whether the 
axino can explain the present DM abundance. There are four production mechanisms of axinos. 
We sketched these four processes in fig]3l One of the processes is the production from the decay 
of the modulus directly. The second comes from the decay of next LSPs (NLSP), Xi produced 
by the modulus. We should note that the sparticles are not in the thermal equilibrium because 
the reheat temperature after the X decay is not high enough. The abundance of axino produced 
by this process is essentially the same as that of NLSP, and hence we estimate the abundance of 
NLSP at the modulus decay. In the third process, axinos are produced by the decay of gravitinos. 
Since gravitinos decay uniformly into all MSSM particles and the axino, the axino abundance 
may be determined by the ratio of the decay channels. The fourth and last process is that NLSPs 
produced by gravitinos decay into axinos. We discuss these processes for cases where the NLSP 
is the bino, the higgsino and the stau in turn. 
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5.1 Bino NLSP case 

1st. process. The branching ratio of X into the axino pair can be estimated from eq. (1511) and 
eq.(l52l as 

fl#W ^ S) .J-^)\ (56 ) 

According to eq.(|56l) and eq.( l53l . the relic abundance of axinos directly produced by the modulus 
decay is 

* - fa* „ 9 , x ™ (^.)- I/4 v (^f. (57 , 

For the axino mass below the electroweak scale, this yield will be too small to account for the 
present DM abundance. 

2nd. process. When the modulus reheats the Universe, neutralinos are so abundant that the 
annihilation among them may become effective. However, since the neutralino can decay into the 
axino, if the decay width is larger than the interaction rate, then the annihilation process will not 
be effective. The bino decays into the axino and a photon or a Z°-boson through the following 
coupling 

aiN 1 - 

where N is the number of the messengers, the field strength of the U(l)y gauge boson B^. 
In this expression, we used the four component notation. Eq. (l58l) leads to 



~ a 2 N 2 1 m 3 ~ 

r(i^a 7 )^ -_ 2Q (59a) 



256tt 3 cos 2 6 w (S) 2 ' 



1 {B aZ j - 256tt 3 cos 2 6 W { 1 XZ) V 2 2 J (5)2 ' ^ 

where a em is the fine structure constant, 6 W the Weinberg angle and %z = m 2 z0 /m~ with Z°-boson 
mass, m z o, and the NLSP mass, m^. When the bino is sufficiently heavy, its lifetime is relatively 
short 



T B 



1.9 x 10" 4 sec., (60) 



for = 300 GeV and (S) = 10 10 GeV 10 . Assuming that the bino is heavy enough to annihilate 
into a top quark pair through s-wave, the annihilation cross section of the bino, ((J axax ,v T ^ l .) §, is 



10 



In the mirage mediation, such a relatively heavy bino is usually realized. 



17 



given by 



327r m\ 




<<Wi*L>5 - ^i (mffl+m |_ m 2 )2 I 1 - -2 1 (61) 

with mt, m t R an d m _g being the masses of the top-quark, the right-handed stop and the bino, 
respectively. Here ffann. is the annihilation cross section of two binos, v Te \, their relative velocity 
and (• • • ) represents the thermal average n . Thus, from eq. (1591) and eq. flST]) . we can find that the 
interaction rate at the modulus decay is much larger than the decay width: 

(Oann.^rcl.) b U B = (^ann.^rel.) § Y B S > ( 62 ) 

where is the number density of the bino. The annihilation process will terminate when the 
Hubble parameter becomes comparable to the annihilation rate 

(Oann.^rel.) B U B ~ H ( T r)- ( 63 ) 

The bino abundance after their annihilation is estimated as 



1 / 90 x 1/2 



T x 4 \ir 2 g*J (cr a xm.Vxel) S T R M pi 



(64) 



After the annihilation, the bino NLSP decays into the axino, and hence the axino abundance is 
the same as that of bino: 



* 6 2 x " /4 d- ■ " n ( mx ) ' m GeV " 2 

T x " V 10 J ' V 12 / \W s GeVJ W m v, d .} s 

' H 9.4X10- (^y'V.ffl^m-' (65) 



yB _ ™B 



10 J 9 V 12 / V10 6 GeV 

where we have set m t = 174 GeV, = 300 GeV and m^/m^ = 1.25 as the reference values. 

3rd. process. A fraction of the gravitinos decays into axinos. Such an axino abundance is 
provided by the ratio of the decay channels. Since the total decay width of the gravitino with 
negligible final state mass is (see, for instance [62J) 

r 244 m 3/ 2 

13/2 = 384^' (66) 



11 If the bino is so light that it cannot annihilate into the top quark, since the annihilation cross section becomes 
proportional to their relative velocity, the cross section is more suppressed. In such a case, the relic abundance 
of the bino becomes more redundant, and hence in order to explain the present DM abundance, the axino mass 
should be smaller about one order or so than that in our case. 
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the axino abundance can be estimated as 

Y~ /2 ~ — Y x , = — -B« - T * 



244 1 3/2 244 2 3/2 



m x 

, , 7xl o-» (B7) 

where Kw 2 denotes the yield of the gravitino produced by the modulus decay In the last step, 
we used the branching ratio into the gravitino pair B 3 / 2 which is defined by 

^ = st®"* a01 - (68) 

4th. process. Finally, let us consider the process that the binos produced by the gravitinos 
decay into axinos. In distinction to the case of the modulus decay, the decay temperature of the 
gravitino is not so high, 

90 N 1/4 



^ 9.0 x IP- 3 GeV f glgVg) V V4 f "°/M 3/a , (69) 
V 10 J Vl0 5 GeV/ ' v ; 

and then, the annihilation process does not occur effectively. Therefore, the axino abundance is 
the same as that of the gravitinos: 



3/2 V 10 / V 12 / V10 6 GeV/ 



(70) 



Comparing eqs.( |57j) . ( |65l) . ( IBTI) and ( 1701 . we find that the axino abundance is dominantly consti- 
tuted from the 4-th. process, that is, decay of the binos produced by the gravitinos, eq. (i701) . 

We should estimate whether the binos can become non-relativistic by scattering with thermal 
bath [63] [61]. At the gravitino decay, binos lose its energy mainly by scattering with the back- 
ground electron. When the bino is produced by the gravitino, it is relativistic, and then the cross 
section, <7 scatt ., is estimated by 



E 2 ~T 2 

(cw^rei.) ^ 128nal^-^, (71) 

e-R B 

where Eg is the energy of the bino and rrie R the mass of the right-handed selectron. The energy 
loss rate for the relativistic NLSP, r^f t p , is given by 

rfcatt. - n e (<Xscatt.'Urel.)^P a , (72) 
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where n e is the number density of the back-ground electron and AE^/Eg the averaged energy 
loss rate of bino in one scattering, which is given by 

AE„ fE n T- in \ , 

^ aU (f} (73) 

Comparing the energy loss rate eq. (172l with the decay width eq.(l59l). we can find 



scatt. 



r 5 



2 x 10 d (74) 



T 3/2 



for rrig ~ rn^ R = 300 GeV, T 3 / 2 — 10 MeV and Eg ~ 50 TeV. Therefore, at the gravitino decay 
the binos decay into axinos after becoming non-relativistic without the annihilation process: 



r s^att. > T B > (<W.«»l) fl "5- (75) 



5.2 Higgsino NLSP case 

Let us next consider the case of the higgsino NLSP. The results of the 1st. and 3rd. processes 
remain valid, irrespective of the composition of NLSPs. A difference from the bino NLSP case is 
the annihilation cross section. It is given by 



arm 



v , )~ ~ 91 1 ( 1 ~^) 3/2 + £ 1 (76) 
reL/?i " 32vrm| (2 - x w ) 2 64tt cos 4 ^ ™f (2 - x z ) 2 1 J 

with xw = wi^/m-, mvi/ being the H^-boson mass, g 2 the SU(2) gauge coupling constant and 
the higgsino mass. Here coannihilation is not taken into account. The higgsino can decay into 
the axino and the Higgs boson through the //-term, eq. (l2"9"|) . 

C = k^4- { /5 )(ah° 1 hi + ah° 2 K) + h.c. 1 . (77) 



The physical neutral Higgs bosons are related to h\ and hi by 

h\ = -j= («i - h° sin a + H° cos a + iA° sin (3 + tG° cos $) (78) 

h° 2 = — (y 2 + h° cos a + H° sin a + iA° cos (3 - iG° sin /3) (79) 
v 2 

with t?i and v 2 being VEVs of the Hi and H 2 . Here H° and A are physical scalars, which 
correspond to the CP even light Higgs, the CP even heavy Higgs and the CP odd Higgs, respec- 
tively, and G° is the Nambu-Goldstone boson. The mixing angle a (which should not be confused 
with a parameter eq.©) satisfies the relation 

tan 2a = ~4 f- tan 2/3, 80 

m A ~ m zo 
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where tan/? = v 2 /vi and is the mass of A . Then, the decay width of the higgsino into the 
axino and the Higgs boson is 

F(h± - ah ) ~ -L (cos a =F sin a) 2 ( 1 - (JL) * ^1 (81) 

32tt V m {J \ m hJ ( S ) 

where \i is the //-term, i.e., eq. (l30l and h± the mass eigenstate, which are defined as h± = 
-^(h\ ± h 2 ). The higgsino can also decay into the axino and the Z° boson through the axino- 
higgsino mixing. According to eq. (l29|) . the axino and the higgsino are mixed each other by 

y2\s 2 / V2\s 2 

This mixing can be removed by transformations 

a' ~a--hS- -hi (83) 

hi' ~hl + -a, (84) 



£ mix = -K—(^J-\v 2 ahl - K^(^-\ Vl ah° 2 + h.c. (82) 



h ' ~h° 2 + -a, (85) 
- n 



I F I \ 

where we have neglected the higher power of ei/ a = — k — ( )Vi/u ~ Vi/(S). In the unitary 

V2\S 2 / 

gauge, the relevant term to the decay is 

£ = -^V-^"^, + ^Ss^flZ, (86) 
2 cos jU ^2 cos W \i 

Thus, the decay width of the higgsino into the axino and the Z° boson is 

T(h± - aZ°) * -L ( C os (3 T sin /?) 2 ( 1 - V ( 1 + 2 ) (87) 



h 

The lifetime of the higgsino can be estimated as 

r~ h = T~ 1 = (r(h -> a/i°) + r(fc -> aZ )) 1 ~ 2.4 x 10~ 9 sec, (88) 

for m h o = 120 GeV, m z o = 91 GeV, m~ h = 150 GeV and (S) = 10 10 GeV. According to eq.flZSJ and 
eq. llHT!) and eq. llHTI) . we can find that in the 2nd. process, the annihilation is more effective than 
the decay. Since the higgsinos decay into axinos after their annihilation, the axino abundance of 
the 2nd. process is comparable to that of the higgsinos 



u Tlr 
yh ^ _h 



t$ \ 10 7 9 V10 6 GeV/ ( 
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Here, in the last step, we have set = 150 GeV for reference. 

On the other hand, in the 4th. process, the annihilation does not occur effectively as in the 
case of the bino NLSP. Therefore, the axino abundance produced by the 4th. process is the same 
as the gravitino abundance 

y} * y 3 % * 4.3 x io- (^ilY 1/l i, (^Y 2 f^M" 2 , (90) 

3/2 V 10 / V 12 / V10 6 GeV/ ' V ; 

and hence we find this is the dominant contribution to the axino relic abundance. As the case of 
the bino NLSP, the higgsinos also become non-relativistic before its decay. The scattering rate 
for the higgsino is given by [61] 



m W V m h m z 



12^-2^|iV f) (91) 



m~ mr 

h n 



where Am^ = — is the mass difference between the chargino and the higgsino, Er the 
energy of the higgsino and Np the number of the processes. Taking the ratio eq. llHTT) to eq. (!9T|) . 
we find 



r 



h 

scatt. 



400A^ F , (92) 



T 3/2 



for E~ h = 50 TeV, m~ h = 150 GeV, Am x ~ ~ 5 GeV and T 3/2 ~ 10 MeV 12 . 



5.3 Stau NLSP case 

Next, we discuss the case for the stau NLSP. The annihilation cross section of the stau is given 
by [65] 

4:iral m 167rag m m| 

(eWUrel./r ^ ~p H —3, (93) 

m r cos 4 9w{m~ + m z ~) 

where is the stau mass. For simplicity, we do not consider the coannihilation of r with B. 
The stau can decay into the axino through the axino-higgsino mixing eq. (|82[) . Thus, the coupling 
of the stau with the axino becomes 

Crra - -Vt—TrT a' + h.C. (94) 



12 In the other NLSP cases, one would also find the scattering process of NLSP with thermal bath is more 
effective than its decay width. Therefore, any NLSPs become non-relativistic before its decay. 
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with y T the Yukawa coupling of r. Therefore, we can obtain the decay width of the stau is 



1 / m T \ m~ 



32vr V^r/ (S) 2 



(S) \ 2 ( TUr Y 1 



The lifetime is estimated as 

T ^2.0xlO-sec.( I ^j . (96) 

In the 2nd. process, since staus are produced by gauginos which are the decay product of the 
modulus, they are so abundant that the annihilation becomes effective. Thus, the axinos produced 
in the 2nd. process can be negligible. Comparing the decay width of the stau with the interaction 
rate by using eq.( l93l) and eq. (195j) . we can find that the annihilation of the stau in the 4th. process 
does not become effective. Thus, when the stau is NLSP the axino abundance is also equal to 
the gravitino abundance, 

5.4 Stop and Wino NLSP cases 

We also mention the case of the stop and the wino NLSP briefly. Let us discuss first the case 
where the stop is the NLSP. If the stop mass is heavy enough to decay into the top quark, the 
decay width of the stop is obtained by (analogous to eq.flS! 



^ If m 2 t \ fmt\ m r 

r(t^at)~ 1 1 — -^S 98 

32ttV ml \m T J 



and hence the lifetime of the stop is 



(3) \ f mj \ 



^ L8xl0 " 8se Hl^J (Solfcvj (99) 

with m t = 174 GeV. The annihilation cross section of the stop is obtained by replacing the 
coupling constant of eq.(l93l with the strong coupling, 

. . 327raq 

(<T ami v Te i.)t^ (100) 

m~ 

Thus, when mj > m t GeV, we can find the annihilation process is not effective in the 4th. process. 
Then, the abundance of the axino is the same as that of the other NLSP cases, i.e. eq.( 170l) 13 . 



13 On the other hand, if mj < mt , the decay width of the stop is more suppressed than eq. ([98|) by the three-body 
phase space. In such a case, since the annihilation may become effective, the axino abundance in the 4th. process 
is much smaller than that of the case. Therefore, when mj < mt the axino cannot explain the present DM 
abundance. 
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If the wino is NLSP 14 , its decay width can be obtained by the analogy of that of the bino: 

a 2 N 2 1 m~ r 
256vr 3 sm 2 6 w (S) 2 



~ a 2 N 2 1 ™%r r 

W - St) - ^zr-^-^ (101a) 



^^l&S&M'-T-^lfa (1 ° lb) 

where is the wino mass. For = 100 GeV, the lifetime of the wino is 

~ 1.9 x 10" 3 sec.. (102) 

On the other hand, it is known that the annihilation of the wino pair is very effective. The 
annihilation cross section is given by [61] 

STraj (1 - x w f' 2 
(a^v Tel ) W = {2 _ xw)2 ■ (103) 

Then we can find that when ~ 100 GeV, (cr ann v Te i)^ a ^ the g ray itino decay, where 

riyy is the number density of the wino. Thus the abundance of the axino is the same as that of 
the wino after annihilation, 



T 3/2 



where we have used xw — 0.64 evaluated by = 100 GeV. On the other hand, when > 
300 GeV, the annihilation process does not occur effectively, and then the axino abundance is 
equal to that of the other NLSP cases. 

Here, we summarize the axino abundance for each NLSP cases in table. 1, where we showed 
only the 4th. process because the other processes are subdominant. When B, h and r are NLSP 
cases, the annihilation is always less effective than the decay. Since the annihilation cross section 
of the stop is quite large, the axino abundance becomes too small to explain the present DM 
abundance when the annihilation of the stop is effective. 



6 Results 

At last, we are now ready to discuss whether the axino relic abundance can explain the present 
DM abundance. If the decay of NLSP is more effective than the annihilation in the 4th. process, 

14 This does not occur for c\ a = 1 and a ~ 1. However, once we introduce e.g. non-universal eg [S7]-[nH], it 
is possible to have wino NLSP, which will not upset the stabilization mechanism of S controlled by the strong 
interaction unless bino is hierarchically heavy. 
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NLSP 


Case A: T ^> n(av) 


Case B: T <C n(av) 


B 






h 




does not occur 


T 


Y~ a ~ 4.3 x 10~ 9 




t 




too small abundance 


w 




Yz ~ 3.3 x KT 11 



Table 1: Table of abundance of the axino in the 4th. process. Contributions of the other processes 
are subdominant. Cases A and B, respectively, correspond to situations where the annihilation 
process of NLSP is not effective and effective. 

the Q parameter, which is defined by the ratio of the mass density of a and the critical density, 
is the same in any NLSP cases 

^ s L2 (ja_) (<±my 1/, dg f^V /2 f^M" 2 . (io5) 

VI GeV/ V 10 / V 12 / V10 6 GeV/ V ; 

On the other hand, if the wino is NLSP and their annihilation is effective, the Q parameter is 

nia \ ( ( g*(T 3 /2) \ ( m 3/2 \ -3/2 / m ~ v 2 



n-ztf ~ 10- d — — W 26 . , , 

.lGeV/l V 10 J V50TeV7 VlOOGeV/ 

x\ \ -1/4 



+ (4 , + , 5A W) (*&!)- (^)-|. (1 oa, 



Here we have set d g — 1 and A^c = 12. In figJU we plotted eq. (jl05j) in terms of mx and 
rrta- From the WMAP three year results [70], the DM abundance in the present Universe is 
Vl DM h 2 = 0.105j£oi£ (68 % C.L.). Thus, we can find from ngHthat the axino with m~ a ~ 0(100) 
MeV can explain the present DM abundance in any NLSP cases, if the decay of NLSP is more 
effective than the annihilation in the 4th. process. In fig El we have plotted axino mass contours 
which satisfy VLah 2 = 0.1 in 7713/2 — mx plane for the wino NLSP case, where we set A = iV = 1, 
k = 0, and = 100 GeV. If the wino is NLSP and their annihilation is effective, figJU shows 
that the right amount of DM can be explained by the axino with a few GeV mass 15 . Here we 
should caution that the estimate of the relic abundance given in this paper is rather rough, which 
may contain an error of factor 2 or so. From eq. (1251) . we can parametrize the axino mass by using 
A as 

m~ a = ^-F*. (107) 



15 Here, we should consider that the gravitino mass is heavier than about 30 TeV not to spoil the success of BBN 
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Substituting eq. (11071) into Q parameters, we can find 



V0.06; VSOTeVy V10 6 GeV/ V 10 J 9 V 12 ' 
for the case where the annihilation is not effective, and 



' O^JUoTevJ 1 l 10 J v50TeW llOOGeW 



for the wino NLSP case, whose annihilation process is effective. Thus, (9(0.1) coupling between 
the axion superfield and messengers leads to the right amount of DM in any NLSP cases. 

Here, we should mention the free-streaming scale of the axino LSP. According to the discussion 
in the previous section, the axino LSPs are produced mainly at the decays of non-relativistic 
NLSPs. At the production, they are relativistic and travel freely until their momentum gets red- 
shifted to be non-relativistic. The observations of Lyman-a forest require that the free-streaming 
scale of DM have to be less than 0(1) Mpc [74] [75]. We can estimate the free-streaming scale of 
the axino as [76] [77] 



A FS = / ^dt~1.0Mpc« a 



a(t) 



10 6 ; 



1 + 0.14 In 



10 6„Nl/2 



ta 1 



(110) 



where v(t) and a(t) are the velocity of the axino and the scale factor, respectively. £eq is the time 
at matter-radiation equality, t% the time when the axino is produced and the three-momentum 
normalized by the axino mass Ua = p/rria- From the discussion of section [51 we can approximate 
that ta is equal to the lifetime of the gravitino, and hence we find 

\ nA\/f ( m NLSP W m~ a m 3/2 \-3/ 2 

for any NLSP cases if the annihilation process is not effective. The axinos produced directly by 
the decays of moduli and gravitinos have much larger free-streaming length of order 10 Mpc. 
However, their contamination in the total dark matter density is about 1% level, so that their 
contribution will be harmless [74J. 

On the other hand, in the case where the annihilation process is effective and NLSP is the wino 
with rriyy = 100 GeV, axinos produced by gravitinos and the modulus have the free-streaming 
scale of the order of 0(1) Mpc 16 : 

16 The free-streaming scale of the axino produced by the stau or the wino NLSP is 0(0.01) Mpc. 
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The free-streaming scales of the axino, eq. fillip and eq. (11 121) . are in an interesting range for the 
small structure problem point of view. The model of the cold dark matter (CDM) can explain the 
large scale structure of the Universe very well. However, iV-body simulations in the CDM model 
seem to conflict with observations at the smaller scale (< Mpc). Numerical simulations based 
on the CDM model predict that the number of halos is by one order of magnitude larger than 
observations within 0(1) Mpc. This is called the missing satellite problem [78]. Another possible 
discrepancy is that the prediction of simulations for mass profile of CDM halos is excessively cuspy 
compared with observations, which is known as the cusp problem |79j . One of the resolutions to 
these problems is to invoke some DM whose free-streaming scale is in the range of 0(0.1 — 1) 
Mpc, because it can sweep out density fluctuations in the small scale [80] . Then the axino LSP 
of our model may also be able to solve these problems. 

Before concluding this paper, we also mention whether the axion can explain DM. We have 
already discussed that the singlet, S, can be stabilized at 10 10 GeV to 10 12 GeV by choosing 
the modular weight of the messenger appropriately. The axion relic abundance eq. fl22|) with 
= 0{1) shows that when the PQ scale is obtained /pq ~ 10 12 GeV the axion can constitute 
DM abundance. When the axion superfield is stabilized at 10 12 GeV, however, the lifetime of the 
saxion is longer than that of the modulus, then the Universe may be dominated by the saxion 
if the initial amplitude of the saxion is not so small compared with the Planck scale. Since the 
saxion mainly decays into axions, it will spoil the success of BBN. Therefore, to explain DM by 
the axion would be difficult in our model unless the initial amplitude of the saxion is suppressed 
or there is an extra entropy production. 

7 Summary 

Let us summarize this paper. We have discussed the axionic extended mirage mediation, axionic 
mirage mediation, to remedy crucial problems included in the mirage mediation: the /i-problem 
and the moduli problem. We showed that the PQ scale can be obtained in the axion window, 
10 9 GeV < / P q < 10 12 ~ 13 GeV, by stabilizing the axion superfield. Especially, in the KKLT 
set-up (a ~ 1) with £ = 1/2, we obtained the PQ scale at 10 10 GeV. To obtain desired values 
of fi- 1 B //-term was achieved by introducing some singlets with couplings eq. (l26l) and eq. fl27|) and 
integrating out the singlet T which becomes heavy due to the VEV of the axion superfield. It 
was also found that the phase of the 5-parameter can be rotated away simultaneously with the 
phases of the gaugino mass and the A-parameter. Thus, SUSY CP problem is also absent in our 
model. 

We investigated the implications of the modulus, whose mass is larger than that of the gravitino 
by several order of magnitude, and the saxion to cosmology. When both the saxion and the 
modulus have the order of the Planck scale initial amplitude, the Universe is dominated by 
the oscillation energy of the modulus eventually. In this case, we found that as long as the 
saxion lifetime is sufficiently shorter than that of the modulus, axions produced by the saxion 
decay would not spoil the success of BBN. Thus, we discussed the cosmic evolution only after 
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the modulus decay. Although there were several axino LSP production processes, the dominant 
contribution to the axino relic abundance came from the decay of the NLSPs produced by the 
gravitinos. We have estimated the axino abundance for the various NLSP cases. However, we 
found that the axino relic abundance is the same in any NLSP cases, if the decay width of NLSP 
into the axino is larger than their annihilation rate. In such a situation, we found that the axino 
mass with 0(100) MeV leads to the present DM abundance. On the other hand, in the case 
where NLSP is the wino and their interaction rate is more effective than their decay process, the 
axino abundance is less than the former case about two orders. In this case, the DM abundance 
can be explained by the axino with ~ (9(1) GeV. When we parametrize the axino mass in 
terms of the Yukawa coupling of S, A, by eq. (1107[) . A ~ O(0. 1) explains the right amount of DM 
abundance, eq. fllOSp and eq. (1109j) . Thus, nevertheless the branching ratio of the modulus to the 
gravitino pair is sizable, our model can solve the moduli problem as well as leads to right amount 
of DM abundance. In addition, we also found that the free-streaming scale of the axino LSP is 
in an interesting region 0(0.1 — 1) Mpc for the small structure problems. 

Finally we would like to briefly mention signatures of neutralino NLSP decays at collider 
experiments. As was discussed in this paper, the axino LSP couples to the Higgs multiplets at 
tree level, whereas its coupling to the gauge sector arises at one loop, and thus is suppressed. 
This implies that the neutralino NLSP dominantly decays to the Higgs boson (as well as the Z° 
boson) as far as it contains a sizable fraction of the higgsino components. We expect that the 
NLSP decay into the Higgs boson will provide a spectacular signal of displaced vertex emitting 
hard jets, if the decay vertex can be reconstructed. Detailed study on the collider signatures will 
be presented elsewhere. 
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Figure 1: The saxion potential in the axionic mirage mediation. The minimal model, / oc X 
[c\ a = 1) is considered with a = 1, k = 0, £ — 1/2, A = 1, and N = 1, which satisfies the mirage 
condition. 



34 



0.5 1 1.5 2 2.5 

a 

Figure 2: The PQ scale /pq ~ (S) in the axionic mirage mediation. The minimal model, f a (xX 
(c\ a = 1) is assumed. Three different cases, I = 0, 1/2, 1 are shown as a function of a for k = 0. 
M mms is the mirage messenger scale at which m 2 s (l = 1/2) crosses zero. The remaining parameters 
are fixed to A = 1 and N = 1, for which / P q shows quite modest dependence. For k > 0, /pq is 
considerably lowered than depicted in small A region (< 1) due to the modulous contribution to 
m%. 
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1st. process 
X > aa 
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Figure 3: The sketch of four production processes of the axino, where \ denotes NLSP. The 4th. 
process gives the dominat contribution to the relice abundance of the axino. 




m x [GeV] 



Figure 4: Contours of the density parameter of the axino, fi^/i 2 , drawn in mx — rria plane. Three 
lines represent flah 2 = 0.2, 0.1 and 0.05, from the above. The contours are the same in any NLSP 
cases if the decay of NLSP produced by the gravitino decay is more effective than the annihilation 
process. 
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Figure 5: Contours of the axino mass, ms, which satisfy VL^h 2 = 0.1 drawn in — mx plane in 
the case where NLSP is the wino with ~ 100 GeV and their annihilation is effective. Three 
lines represent = 1, 3 and 5 GeV, from left to right. Here, we set A = N — 1, k — and 
= 100 GeV for simplicity. From the gravitino problem point of view, the gravitino mass should 
be heavier than about 30 TeV. 
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